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MSTC  Organization  &  Activity 


Mission:  Integrate  multiple  disciplines  to  discover  and  exploit  new  phenomena 
for  system  optimization  and  assessment  of  revolutionary  aerospace  vehicles 


Center  Director 


Branch  Chief 
-  Tech  Advisor 


Prototype  Representation  & 

Analysis  Methods  for 

Prototype  Validation  & 

Design  Exploration  Methods 

Prototypes 

Assessment 

•  Parametric  Geometry  & 

•  Multidisciplinary  Analysis 

•HiFi  QTA 

Mesh 

•Appropriate-fidelity  Solutions 

•  Prototype  Experimental 

•  Subsystem  Representation 

and  Sensitivities 

Validation 

•  Design  Space  Exploration  & 

•  Nondeterministic  Models 

•  TRL  Assessment 

Optimization 

•  Risk-based  Design 

Shared  Activity  -  Utilize  a  Unified  Framework  (SORCER,  MODEL  Center) 

-W 
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Some  Significant  Collaborations 


MSTC  Collaborative  Center  with 
VPI  &  SU,  WSU,  and  University  of  Maryland 
(Formed  March  2009) 

Prof  Kapania,  Director 
Dr.  Kolonay,  PM 


WRIGHT  STATE 


VirginiaTech 


Invent  the  Future 


AFRL/RB  and  WSU  Center  for  Micro  Air 
Vehicle  Studies  (Formed  June  2010) 


Prof  Missoum,  Mr.  Basudhar  (UA,  Tucson)  and  Dr.  Lambe  (MSSRC)  -  RBDO  with  LCO 
Prof  Dong  and  Mr.  Gaston  (WSU)  -  ROM  and  Simulation  of  falling  bodies 
Prof  McFarland  and  Mr.  Hubbard  (UIUC)  -  Transmission  design  with  nonlinearity 
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Internal  Collaborations  in  MAVs 


Math  (6.1)  Physics  (6.1) 


Risk-Based 

Computational 

Prototyping 

Beran  (PI), 
Camberos,  Lindsley 


Physics-Based 
Design  Analysis  of 
MAVs 

Snyder  (PI),  Beran, 
Kolonay 


NRC:  Chabalko,  Kurdi,  McClung,  Stanford 


V. 


Basic  Research  in  Computational  Design 
(2009-2011) 

- v - 


J 


Flapping  Sciences  Integration 
(2009-2011):  6.2 


•  Service-oriented  framework 

•  In-house  computer  scientists 

•  Design  tools  (Transition) 

•  Funded  follow-on  design 
program  (MPP,  FY12+) 


AFRL/RB 

MAVTeam 


Basic  and  Applied 
Research  in  MAVs 

•  Structural  and  Flight 
Testing  (Parker)  - 
validation  of  structural 
and  system  models 

•  CFD  (Visbal)  -  verification 
of  aero  models 

•  Controls  Science  (Doman) 
-  integration  of  controls 
models 

•  Unsteady  Aerodynamics 
(OL)  -  validation  of  aero 
models 

•  Perching  Technologies 
(Reich)  -  application  of 
aero  models 
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Role  of  Computational  Mathematics 


PF  =  Probability  of 
Failure  (large  LCO) 
Must  be  sufficiently 
small 


PDF 


Distribution  of  critical  speeds  arises  via 
variability  in  air  vehicle  or  environment 


Large  Vehicles 


Speed  at  which  LCO 
*  achieves  maximum 
value 

►  u 


oo 


,  Subcritical 


Acceleration 
Limit  — . 

/ 

T  j 

Strong  LCO 

Risk  Region  \ 

1  i 

— : — \ 

P - *  Mach 

Flutter  Point  Limit 
(Linearized  Analysis) 


Computational  mathematics  needed  for  physics-based  design  of  reliable  vehicles 
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Role  of  Computational  Mathematics  (cont.) 


•  Exploit  nonlinear  aeroelastic  interactions  for  small  aircraft 


Unsteady  Flow,  lida  (2004)  Unsteady  Deformations 


•  Numerous  challenges  for  design  of  Micro  air  vehicles  (MAVs) 

-  Physics  Rich  (must  be  a  physics-based  approach) 

-  Complex  and  time-dependent  actuations  (unsteady) 

-  Non-conventional  geometries  and  structural  topologies 

-  Power-based  integration  of  propulsion,  structure,  control  components 

Computational  mathematics  needed  for  physics-based  design  of  MAVs 

Akiyoshi  lida  (2004)  -  Kogakuiu  Univ,  http://www.nagare.or.jp/mm/2004/gallery/iida/dragonfly.html 
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Spectral  Formulation  for  Time-Periodic  Systems 


Requires 

increased 

resolution 


Monolithic-time  collocation 


Attack  high-frequency  content 
with  ftp-refinement 


Attack  Large  #  of  space  DOFs 
with  ROM 


Dimensionality  of  X  X 


Number  of  DOFs  in  time 
(elements  x  order+1) 


X, 


mon 
Total  DOFs 


Uses  a  local  basis  instead  of  global  basis 

m 

x,(0  =  2x.((,)'r«(0 

q= 0 

M  -  Order  of  the  spectral  element 

£q  -  Zeroes  of  the  Lobatto-Legendre  polynomials 

xj/  (g)  -  Lagrange  polynomial  of  order  m 


Kurdi  and  Beran,  "Spectral  Element  Method  in  Time  for  Rapidly  Actuated  Systems JCP,  Vol.  227,  No.  3,  2008,  pp.  1809-1835. 

7 

Approved  for  public  release. 


Monolithic-Time  Collocation 


Arrays  corresponding  to  a 
discrete  2D  field  variable 


X  t°  t1  t2  t3 


Xmo„  =  [X°,  X1,  X2,  X3, ...] 


Context  for  time-periodic  and  transient  solutions 
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Adjoint-Variable  Approach 


Solve  Fmon(Xmon,X)  =  0 


adjoint 


H(Xmon)  =  objective 
Fmon  =  equation  residual 


Sensitivity: 


dH 


dH 


dx  ax 


0F. 


\ 


-1 


mon 


ax. 


ap 


mon 


mon  \  u  'mon  / 

v _ 


dX 


J 


direct 


O 


O 
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Adjoint  Computation  for  Transient  Sensitivity  Analysis 


Goal:  Develop  a  sensitivity  analysis  process  that  scales  well  with  total  #  DOFs 


•  Interested  in  the  adjoint-variable  approach  in  anticipation  of: 

-  many  design  variables  (not  true  of  direct  and  sampling  based  approaches) 

-  use  of  gradient-based  optimization  (trade  global  effectiveness  for  efficiency) 

•  Some  relevant  literature 

-  Nadarajah  and  Jameson,  "Optimum  Shape  Design  for  Unsteady  Flows  with 
Time-Accurate  Continuous  and  Discrete  Adjoint  Methods,"  AIAA  Journal  Vol. 
45,  No.  7,  2007 

-  Thomas,  Hall,  and  Dowell,  "A  Discrete  Adjoint  Approach  for  Modeling 
Unsteady  Aerodynamic  Design  Sensitivities,"  AIAA  2003-0041,  2003 

-  Mani  and  Mavriplis,  "An  Unsteady  Discrete  Adjoint  Formulation  for  Two- 
Dimensional  Flow  Problems  with  Deforming  Meshes,"  AIAA  2007-60,  2007 

•  Create  a  sample  problem  to  explore  a  POD-based  approach  to 
eliminate  challenge  of  storing  the  forward  solution 
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Problem  Description 


Transient  analysis  of  incompressible  flow  in  a  square  cavity  with  unsteady  lid 


y 


Flow  at  rest  initially 


•  Steady:  U  =  1  (impulsive) 

-  verify;  assess  accuracy 

•  Transient:  U  =  >z(l-cos(f  t)) 

-  define  H,  a  function  of  the  transient  solution 

-  compute  sensitivity  of  H  to  frequency,  f 

•  Streamfunction-vorticity  form 


dco  dco  dco  1  2  o  dV 

—  +  u —  +  v —  =  — V  to  V  ^  =  -00  u  =  — -  v  =  — — 
dt  dx  dy  Re  dy  dx 


11 

Approved  for  public  release. 


Discretization  and  Time  Integration 


Explicit/implicit  formulation 


2nd-order-accurate,  centra  I -difference  operators 


Repeat  for  next 
time  step 


co(x,l) 


4('I,(x,l-Ay)  +  U(t)Ay)  +  0(Ay) 
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Adjoint-Variable  Approach 


Linear,  time  invariant 


n  =  1 


n  =  2 


n  =  3 


0 

0 

0 

0 

"CO 

I 


0 

0 

0 

0 

0 

L 


V 


^^mon  ^mon 


Reverse-time 

A 

r| 

0 


Vorticity  BC  coupling 
terms 


Jacobians  arising  from  convective 
terms  [apply  data  compression] 


0 

0 

0 

0 


0 

I 

0 

0 

0 

0 


dH 


v 


ax 


mon 


J 


Linearization 


Verification  (Steady  State) 


tfjna|  =  20,  time  step  =  0.001  tfjna|  =  50,  time  step  =  0.001 


Re  =  100 


Re  =  1000 


*Sahin  and  Owens,  "A  Novel  Fully  Implicit  Finite  Volume  Method  Applied  to  the  Driven  Cavity  Problem  -  Part  I: 
High  Reynolds  Number  Flow  Calculations/'  IntJNum  Methods  Fluids,  Vol.  42,  Issue  1,  May  2003,  pp.  79-88 
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Verification  (Transient) 


Current 


Sahin  and  Owens  (2003) 


•  Re=10000 

•  U(t)  =  1 

•  Contour  plots  of  XV 

•  t=2:  agree  within  2.8% 

•  t=8:  agree  within  4.4% 

•  Need  to  explore  mesh 
and  time  step 
refinements 
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Verification  (Sensitivity) 


t=2  t=4 


U=0.71  U=0.83 


t=10 


U=0.92 


•  Re=1000  with  baseline  mesh  (101  x  101) 

•  U  varies  in  time 

•  Determine  sensitivity  of  H2  about  f  =  1 

•  H2  evaluated  at  t  =  10 

•  Finite-difference  sensitivity:  6f  =  0.0001 

•  Sensitivities  match  to  6  significant  digits 


U(t)  =  l(l-cos(ft)) 

H2(Xmon)  =  2K)2 


dH2/df  (Adjoint) 

dh2/df  (Finite  Difference) 

4.70771958780 

4.7077182309 
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POD  Data  Compression  for  Sensitivity  Analysis 


•  Same  conditions  as  verification  case 

•  Integration  time  of  10;  1000  time  steps 

•  Collect  snapshots  once  every  10  time  steps 

•  Decimate  snapshot  set  to  coarsen 

•  Evaluate  efficiency  and  accuracy  of  POD-based 
adjoint  sensitivity  analysis  as  function  of 
number  of  snapshots  and  modes 


t  =  0  to  10 


t  =  0  to  10 


Cavity 

snapshots 

Compute 

modes 

Cavity 

analysis 

■  W\ 

POD  modes 

■  W\ 

analysis 

On-the-fly,  use  the 
modes  and  the 
instantaneous  modal 
amplitudes  to 
reconstruct  the  flow 
solution 


modes  modal 

amplitudes 

t  =  10  to  0 


Linearized,  adjoint-variable 
sensitivity  analysis 


m 
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Solution  and  POD  Modes  (Streamfunction) 


Solution 


Mode  1 


Mode  2 


■ -0.002 
I-qcw 
1-a.me 
1-0.006 
H-ooi 

I  -O.O'fS 
I  -0  0t4 

1  -aote 
H  -aois 
H-0.02 
|j  -0.022 

g-0D24 

1-0.026 
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Efficiency  and  Accuracy 


dWjdi  using  100  snapshots 


Full  order 

50  modes 

20  modes 

10  modes 

5  modes 

4.707719587 

4.707725353 

4.711403732 

4.724862963 

3.121007234 

%  Error  in  dH2/df 


50  modes 

20  modes 

10  modes 

5  modes 

100  snapshots 

0.00012 

0.078 

0.36 

-34 

20  snapshots 

— 

1.4 

3.2 

-31 

10  snapshots 

— 

— 

1.6 

2.8 

20  snapshots  =  2%  of  time-history  data 
10  modes  =  1%  of  time-history  data 


High  efficiency 
Good  accuracy 


Greatly  decrease  memory  requirement  at  2x  cost:  explore  other  POD  uses 
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Structural  Design  (Inertial  Loads  Only) 


Goal:  study  transient  sensitivity  analysis  in  context  of  DOF  reduction 


Extension:  Kinematic/Structural  design  for  wing  design 
with  airloads  (Stanford  et  al.,  SDM10) 


baseline  - 

optimal"- 

rigid- 


to  =  7Hz 


Identify  best  thickness  distribution  for  rapidly  actuated  plate 

Nonlinear  modeling  of  a  flapping  plate 

Minimize  6ave  =  time-averaged  6 

256  variables  (element  thicknesses  w/  constraints) 

GBO  via  MATLAB  (fmincon) 


Resultant  thickness  distribution 


0,004 


0.003 


Time  Integration  (TI) 
Reduced  Order  (TI) 
Spectral  Element 


5ave/L  0 


002 


5X  design 
speed-up 


0,001  - 


Normalized  design  cost 


ROM  Adjoints  ~  Free 


10 


Stanford,  Beran,  and  Kurdi,  "Adjoint  Sensitivities  of  Time-Periodic  Nonlinear  Structural  Dynamics  via  Model  Reduction/' 
Computers  and  Structures  (to  appear),  2010. 
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lip  amplitude 


Beam  Design  (Inertial  Loads  Only) 


Identify  best  area  distributions  for  minimum  and  maximum  time-averaged  tip  displaced 
Co-rotational  FEA  formulation;  50  beam  elements,  each  with  a  different  sectional  area 
Side  constraints  on  area;  GBO  via  MATLAB  (fmincon) 

Compute  sensitivities  with  the  adjoint  formulation 


min. 


0.4 


0.35 


0.3 
0.25 
60  0.2 


0.15 


0,1 


baseline 


0.2 


0.4  0.6 

x/L 


0.8 


— e —  Newmark,  FOM 
— 6 —  Newmark,  ROM 

— 0 —  SE,  FOM  — e —  FD,  FOM 
SE,R0M  — * —  FD,  ROM 

Some  benefits  of 
combining  SE  and  ROM  in 
time-periodic  formulation 
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Reliability-Based  Design  Optimization  (RBDO) 


Goal:  Examine  use  of  transient  sensitivity  analysis  to  design  a  plate  wing  that  is 

both  light  and  reliable 


Minimize  mass  of  plate;  constrain  the  probability  that  LCOamp  >  6  (PF  <  a) 
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Contrasting  Approaches 


Deterministic  Optimization 


Generally,  the  designed  plate  "moves"  to  the  constraint 


min  weight  =  f(d) 
d 


g(x(d,  E,  IVU)>0; 
side  constraints  on  d 


Generally,  the  designed  plate  "moves"  away  from  the 
constraint  boundary  a  "safe"  distance  (PF  =  a) 


min  weight  =  f(d) 
d 

subject  to: 

1  -  Prob( g  <  0  )/a  >  0;  side  constraints  on  d 


Allen  and  Maute,  "Reliability-based  design  optimization  of  aeroelastic  structures,"  Structural  and  Multidisciplinary 
Optimization,  Vol.  27,  2004,  pp.  228-242.  (Static  Aeroelasticity) 
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RBDO  Formulation 


•  IVL  and  E  are  chosen  to  be  uncertain  (normal) 

•  Map  to  uncorrelated  random  variables  u1  and  u2 
in  standard  normal  space 

•  Compute  Most  Probable  Point  (MPP)  and 
reliability  index  P 

•  Approximate  failure  surface  as  linear:  First  Order 
Reliability  Method  (FORM) 

•  Compute  probability  of  failure,  PF  =  PF((3) 

•  Meet  PF constraint  using  analytical  gradients 


1  For  a  given  structure,  compute  MPP 
using  gradient  based  optimization: 
require  sensitivities  of  g  to  u1  and  u2 

2  Reduce  weight  while  meeting  PF 
constraint  using  gradient  based 
optimization:  require  sensitivities  of 
PFto  d,  (found  from  sensitivities  of  g 
to  d,) 


Adjoints  of 
transient 
solutions  used  to 
compute 
sensitivities  of  g 
to  d, 
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Mach  Range 


RBDO  and  SVM  Results 


Uniform  (baseline)  panel 


M. 


1.  Basudhar  used  Support  Vector 
Machine  and  adaptive  sampling 
to  approximately  construct 
failure  surface 

2.  Computed  PF  with  MCS  on  SVM 
boundary  (55  samples) 

3.  Computed  PF  with  QMCS 
(Lambe,  MSSRC) 

Method  PF 

FORM  0.0197 

MCS  (106)  0.0248 

QMCS  (104)  0.0244 


baseline:  mass  =  0.252  kg 
LCO  =  5.98  mm,  Pf=  0,114% 

amp 


Basudhar  and  Missoum,  "Update  of  explicit  limit  state  functions  constructed 
using  Support  Vector  Machines"  AIAA  2007-1872,  April  2007. 


121 


m 


21\l 


25 


-H 


Mill  i08 


\ 


1,2 


deterministic  optimum:  mass  =  0. 1 78  kg 
LCO  =10.01  mm,  Pf  =  50.17% 

amp 


RBDO  Step  Cost  (MATLAB,  single  CPU) 

Simulation  10  minutes 

Adjoint  5  minutes 

MPP  1  hour 

Optimization  4  hours  (deterministic),  12  hours  (probabilistic) 


probabalistic  optimum:  mass  =  0.182  kg 
LCO  =  7.17  mm,  Pf=  0.130% 

forward 
mass 
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Recent  Activities:  Rigid-Body  MAV  Motions 


Start  to  investigate  impact  of  rigid-body  motion  on  MAV  performance 
Prof.  Haibo  Dong  (WSU),  Mr.  Zachary  Gaston  (WSU) 

Mr.  Tim  Broering  (UL) 


OILS 


UNIVERSITY  OF 

LOUISVILLE 


WRJGHT  STATE 


Chakravarthy,  Albertani,  Evers,  "In-Flight  Dynamically  Adaptive  Configurations:  Lessons  from  Live  Lepidoptera," 
AIAA  2010-2828,  April  2010. 


Digital  Image  Correlation  by  Prof. 
Albertani  using  live  specimens  of 
Lepidoptera 


Camera  1  Camera  2 


McGuire  Center  for  Lepidoptera  and 
Biodiversity,  Gainesville,  FL 


\  /h 


|  ADdomen  root  j 


i  uJV  \  l\i\  n  n  m  m 

tix  WWW  UunMA 


A 


\  AjfvA/vvV 


.  X  A  /I  , 
W\Jw\  /  M 


Need  to  include  rigid-body  motions  and  body  flexibility  in  bio-inspired  MAV  models 
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Plan  for  Rigid-Body  Coupling 


•  Emphasize  passive  motions  first:  falling  bodies  in  quiescent  flow 

-  Pesavento  and  Wang,  "Falling  Paper:  Navier-Stokes  Solutions,  Model  of 
Fluid  Forces,  and  Center  of  Mass  Elevation",  PRL,  Vol.  93,  No.  14,  2004 

-  Modify  high-fidelity  tools  to  repeat  2D  simulations  and  extend  in  3D; 
validate  at  WSU  with  high-speed  photography  (want  comparisons) 

-  Calibrate  quasi-steady  models  (like  those  used  in  flapping) 

•  Examine  influences  of  gust  and  variability  on  falling  motions 

-  Introduce  variability  into  quasi-steady  models  (e.g.,  how  is  seed  dispersal 
impacted  by  winds?) 

•  Re-examine  design  procedures  that  have  been  developed  so  far: 
want  MAVs  that  are  robust  to  gust 
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Some  Typical  Motions 


c  -  chord 
h  =  thickness 
Qf=  water 

\  P 

ps  =  aluminum  s 

ellipse 


Pf 

A 

-gJ 


•  Will  explore  impact  of  variability  (physical 
and  model)  on  distribution  of  landing 
locations 

•  Developed  transient  sensitivities  (direct): 
role  in  selecting  bodies  with  more 
desirable  falling  characteristics? 


tumble 
w/  wind 

flutter 
no  wind 

u. 


7 - 

flutter 
w/  wind 


flutter:  (3  =  h/c  =  1/14 
tumble:  (3  -  h/c  =  1/5 


wind 

10  cm/sec 

- > 


Trajectories  computed  with 
quasi-steady  aerodynamics 
(Wang  et  al.) 


tumble 
no  wind 

(descent  speed  = 
11.5  cm/sec) 


28 

Approved  for  public  release. 


High-Fidelity  Results 


Re  =  40  (Stationary) 

Overture 


N  UNIVERSITY  OF 

LOUISVILLE 


Re  =  01103]  (Falling) 


Preliminary  VICAR3D  result 


WRIGHT  STATE 
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Recent  Activities  (cont.) 


McFarland  and  Hubbard 

Start  to  explore  role  of  actuation  mechanism  in  MAV  design 


-  Investigate  physical  interactions  between  a  flapping  wing  and  the 
mechanism  that  flaps  the  wing  (e.g.,  transmission  of  inertial  loads) 


*1 

mi 

*2 

m 

k 

/  l 

/  | 

/  =0 

/  8 

/  *E 

/ 

..  -I 

Cl 

C'i 

c 

w  -1 

k  r  t  st 

It*  —  —  _  _  _  j 

r  —  “ -  — i 

—►«,•(<)  | - ►  »(/)  | - *>2(0 


•  Developing  compliant  mechanisms  via  topological  optimization 

-  Link  mechanism  with  generated  inertial/aero  loads  (MAO  2010) 


Understanding/modeling  energy  transfers  between  mechanism  and  wing  critical 
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Concluding  Remarks 


•  Sensitivity  analysis  of  transient/time-periodic  systems  serves  an 
important  role  for  design  of  both  large  and  small  aircraft 

-  Constraint  boundaries  often  nonlinear  (LCO  and  aeroelastic  response  in 
gust);  strive  for  physics-based  approaches  not  reliant  on  safety  factors 

-  Essential  for  design  of  flapping  wing  MAVs;  strive  for  physics-based 
approaches  that  account  for  gust 

•  Lessons  learned  through  unsteady  sample  problems 

-  POD  is  a  straightforward  means  for  data  compression  in  sensitivity 
analysis  for  large  systems;  extensions  using  POD  ripe  for  study 

-  Adjoint  vectors  in  ROM  formulation  computed  virtually  for  free 
(tailoring  of  structure  for  nonlinear  response  during  rotary  actuation) 

-  Adjoint-based  sensitivities  work  well  in  an  RBDO  context;  want  to 
extend  (e.g.,  transonic,  SVM,  SORM)  based  on  lessons  learned 

•  Interesting  departure  points  for  further  study:  variability  in 
motion  subject  to  gust,  mechanism  design 
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Recent  Publications 


•  Stanford,  B.,  and  Beran,  P.,  "Adjoint  Sensitivities  of  Time-Periodic  Nonlinear  Structural 
Dynamics  via  Model  Reduction"  Computers  and  Structures  (to  appear),  2010. 

•  Stanford,  B.,  and  Beran,  P.,  "Analytical  Sensitivity  Analysis  of  an  Unsteady  Vortex  Lattice 
Method  for  Flapping  Wing  Optimization,"  Journal  of  Aircraft,  Vol.  47,  No.  2,  Mar. -Apr.  2010,  pp. 
647-662. 

•  Ghommem,  M.,  Hajj,  M.R.,  Pettit,  C.L.,  and  Beran,  P.S.,  "Stochastic  Modeling  of  Incident  Gust 
Effects  on  Aerodynamic  Lift,"  Journal  of  Aircraft  (to  appear),  2010. 

•  Missoum,  S.,  Dribusch,  C.,  and  Beran,  P.,  "Reliability-Based  Design  Optimization  of  Nonlinear 
Aeroelasticity  Problems,"  Journal  of  Aircraft,  Vol.  47,  No.  3,  May-June,  2010,  pp.  992-998. 

•  Kurdi,  M.,  Beran,  P.,  Stanford,  B.,  and  Snyder  R.,  "Optimal  actuation  of  nonlinear  resonant 
systems,"  Structural  and  Multidisciplinary  Optimization,  Vol.  41,  No.  1,  Feb.  2010,  pp  65-86. 

•  Pettit,  C.L.,  Hajj,  M.R.,  and  Beran,  P.S.,  "A  Stochastic  Approach  for  Modeling  Incident  Gust 
Effects  on  Flow  Quantities,"  Probabilistic  Engineering  Mechanics,  Vol.  25,  Issue  1,  Jan.  2010, 
pp.  153-162. 

•  Stanford,  B.,  Beran,  P.,  and  Kurdi,  M.,  "Model  Reduction  Strategies  for  Nonlinear  Beams 
Subjected  to  Large  Rotary  Actuations ,”  Aeronautical  Journal,  Vol.  113,  No.  1150,  Dec.  2009,  pp. 
751-762. 
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Questions? 
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AMP  Team  Composition  (WPAFB) 


Mission:  Integrate  multiple  disciplines  to  discover  and  exploit  new  phenomena 
for  system  optimization  and  assessment  of  revolutionary  aerospace  vehicles 


Center  Director 


Branch  Chief 
-  Tech  Advisor 


Analysis  Methods  for  Prototypes 

•  Dr.  Jose  Camberos  -  On  detail  as  RB  Deputy  Chief  Scientist 

•  Dr  Chris  Chobalko  -  Postdoc  (NRC,  UTC) 

•  Dr.  Ned  Lindsley- Supporting  prototype  validation/assessment 

•  Dr  Aaron  McClung  -  Civil  Servant ,  formerly  NRC 

•  Mr.  John  Moore  -  Undergraduate  Co-op  (University  Florida) 

•  Mr.  Michael  Robbeloth  -  Computer  Scientist,  DSA 

•  Dr.  Rich  Snyder 

•  Dr  Bret  Stanford  -  Postdoc  (NRC) 

•  Dr.  Phil  Beran  -  Lead 
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Methods  Development  Strategy 


Develop  methods:  start 
with  low-dimensional 
formulations  and  move 
towards  high-dimensional 


Goal:  Multifidelity  framework 
built  on  new  methods 


Development  of 
New  Methods 

•  Time-Periodic 
Analysis 

•  Sensitivity  Analysis 

•  Reduced  Order 
Modeling 

•  Uncertainty 
Characterization 


Application  and 
Extension  of  High- 
Fidelity  Methods 

•  Navier-Stokes 
(OVERFLOW) 

•  Beams,  Plates, 
Shells  models 

•  Aeroelasticity 

•  Vortex  methods 
(medium  fidelity) 


Validation 
through  physical 
experiment 

•  Water  channel 
(OL,  AVT-149) 

•  Free  flight  (TU 
Delft,  AVT-184) 

•  Aeroelastic 
ground-test 
facility  (Parker) 


Characterize  physical  Assess  validity  of  all 

limitations  of  lower  methods 

fidelity  approaches 
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Application  to  Insect  Wing 


Berman  and  Wang,  "Energy-Minimizing  Kinematics  in  Hovering  Insect  Flight "  JFM,  Vol. 
582,  2007  (Rigid  wing  with  stroke-plane  deviations) 


Feathering  angle 


Power  reduction  from  initial  design: 

•  55%  for  unconstrained  acceleration 

•  40%  for  constrained  acceleration 
Looking  at  inertial  power  contribution 


Kurdi,  Beran,  Stanford,  and  Snyder,  "Optimal  Actuation  of  Nonlinear 
Resonant  Systems,"  Structural  and  Multidisciplinary  Optimization, 
Published  online  June  2009. 


Prescribed  (cp)  and 
realized  (n)  angles 

•  mass-spring-damper 

•  inertial  &  aero  loads 


Large  snap  rotations  favored 

I  ,\ 

/  X 


t  1;  Optimized  fruitfly 

W  wing  kinematics 

° 7 _ _  (235  Hz) 

-0.2  0  0.2  0.4  0.6  0.8  1  1.2 

Time 
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High-Fidelity  Analysis 


Planform 

Study 

Fx  (N) 

Fy  (N) 

Fz(N) 

Rectangular 

Current  Work 

7.24e-04 

-1 ,46e-03 

7.26e-03 

Manduca  sexta 

Current  Work 

8.42e-04 

-1 ,65e-03 

6.16e-03 

Agrius  convolvuli 

Aono  and  Liu  [1] 

1.20e-03 

-1.20e-03 

8.48e-03 

Understanding 
Complex  Physics 

•  Study  Hawkmoth  physics  using 
Navier-Stokes  (NS)  simulation 

•  Collaboration  with  AFIT 

•  Hawkmoth  kinematics  (hover) 

•  What's  new? 

■  OVERFLOW  2.1  Elastic 
(5th/2nd-order  in 
space/time) 

■  Prescribed  wing 
deformations 

■  Variations  in  kinematics 

•  Moderate  flexibility  increases 
hover  efficiency 


McClung,  "Influence  of  Structural  Flexibility  on  Flapping  Wing  Propulsion,"  AFIT  Dissertation,  April  2009 
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